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DISTRIBUTED  BOSONIC  STATES  AND  CONDENSED 
MATTER  FUSION 


INTRODUCTION 

This  paper  presents  a  theory  of  condensed  matter  fusion1  in  crystals  based 
on  the  periodic  potential  provided  by  an  ordered  lattice  and  the  approximate 
Bose  statistics  obeyed  by  bosons  on  length  scales  associated  with  electrostatic 
interactions.  The  resulting  nuclear  interaction  proceeds  in  a  manner  that 
approximately  preserves  either  the  initial  periodicity  of  the  lattice  or  the 
periodicity  of  a  commensurate  lattice.  Reaction  reversibility  over  timescales 
associated  with  electrostatic  interaction  results  in  condensed  matter  fusion 
being  dominated  by  interactions  in  which  both  the  products  and  reactants  obey 
Bose-Einstein  statistics.  As  a  result  tritium  release  and  3He++  production 
should  not  appear  as  primary  fusion  products. 

BOSE  BLOCH  CONDENSATE 

The  key  attribute  of  Bose  particles  needed  is  the  property  that  the  many 
body  boson  wave  function  possesses  positive  parity  with  respect  to  exchange  of 
single  particle  wave  functions.  As  a  result  any  number  of  bosons  possessing 
a  given  energy  can  occupy  a  particular  lattice  site  (which  is  not  the  case  for 
fermions).  Huge  degeneracy  can  occur  within  the  many  body  boson  wave  function. 
A  condition  for  fusion  is  that  a  portion  of  the  D  atoms  (and  possibly  other 
atoms  and  ions)  that  enter  the  electrode  reside  in  a  screened  ion  form  and  see 
a  periodic  potential.  Such  ions  can  form  a  new  bosonic  state,  which  we  call 
a  boson  Bloch  condensate  (BBC).  The  BBC  is  a  result  of  the  lattice  dependent 
interplay  of  the  symmetries  associated  with  periodicity  and  the  huge  degeneracy 
of  the  boson  wavefunction. 

In  accord  with  Bloch's  theorem2,  the  eigenstates  of  the  crystal  potential 
can  be  written  in  terms  of  single  partied  wave  functions,  which  have  an 
associated  wave- vector  &,  such  that 

iMr+RJ«  exp(iJs  RJ  ^(r)  .  (1) 
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Index  1  designates  the  eigenvalue  ^  associated  with  wave  vector  k  (i.e., 

( k)  «  the  i**  "bosonic  energy  band"  evaluated  with  wave  vector  k).  R„  is  a 
Bravais  lattice  vector.  The  equation  need  not  be  restricted  to  three 
dimensions:  for  the  desired  overlap  to  occur,  the  periodicity  of  the  lattice 
may  be  one-,  two-,  or  three-  dimensional.  Equation  (1)  leads  to  the  obvious 
relationship, 

*i(r)  *i(r+R»)  -  exp(ik  Rn)*V>J(r)  V>i(r)  ,  (2) 
where  tfj(r)  and  ^(r)  are  separate  initial  state  single  particle  wave  functions. 

Equation  (2),  though  seemingly  trivial,  has  important  physical 
consequences.  The  equation  implies  that  even  at  macroscopically  large 
distances  the  overlap  between  boson  single  particle  wave  functions  is  order 
unity.  This  overlap  is  astronomically  large  in  comparison  with  the  overlap 
between  charged  bosons  in  free-space,  where  even  at  molecular  length  scales, 
the  exponential  damping  associated  with  negative  values  of  the  kinetic  energy 
leads  to  entirely  negligible  wavefunction  value.  Because  the  left  side  of  Eq. 
(2)  is  one  of  many  equivalent  single  site  contributions  to  the  many-body 
initial  state  wave  function  (which  involves  comparable  products  throughout  the 
lattice),  it  follows  that  the  fusion  rate  for  bosons  within  the  lattice  should 
be  very  different  than  in  free  space.  This  fusion  rate  is  proportional  to  the 
square  of  a  matrix  element  derived  from  the  integral  of  a  nuclear  potential 
(associated  with  the  strong  forces  that  bind  nucleons  inside  the  nucleus) 
a V,tron*  multiplied  by  the  product  of  suitable  initial  and  final  states.  Because 
the  bosons  are  indistinguishable,  they  behave  as  a  coherent  entity  and  a 
remarkable  effect  becomes  possible:  when  each  ion  responds  to  changes  in  its 
own  local  potential,  it  automatically  induces  changes  in  ill  of  the  potentials 
associated  with  the  remaining  indistinguishable  ions.  Thus,  trough  "self¬ 
interaction"  (the  process  in  which  an  ion  sees  its  own  potential  and  readjusts 
its  motion  and  spatial  distribution,  thereby  modifying  its  own  potential)  all 
of  the  ions  collectively  interact  with  each  other.  Then,  under  suitable 
conditions  (as  discussed  below),  they  can  participate  collectively  in  nuclear 
interactions  at  isolated  locations,  circumventing  the  formidable  tunnelling 
problem  of  free  space  nuclear  physics. 
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Because  the  BBC  is  composed  of  bosons,  eA  may  equal  ei  in  Eq.  (2).  The 
resulting  possibility  of  huge  degeneracies  associated  with  the  many  body  boson 
wave  function  has  important  consequences.  In  particular,  formation  of  a  D+  BBC 
is  a  requirement  for  condensed  matter  fusion.  Such  a  condensate  resembles  a 
quantum  fluid,  such  as  Hell,  in  being  composed  of  indistinguishable  particles. 
However,  the  BBC  is  a  distinct  matter  state,  differing  in  important  ways  from 
such  quantum  fluids.  Unlike  Hell,  BBC  requires  a  host  lattice.  Also,  it  can 
exist  at  very  low  boson  concentration  g  (provided  the  number  of  bosons  NB  is 
macroscopic)  and  does  not  form  exclusively  through  the  ground  state  but 
possesses  many  forms  (or  phases)  associated  with  the  occupation  of  different 
quasi -continuous  energy  bands.  Thus,  it  may  exist  at  elevated  temperatures. 
Also,  the  BBC  couples  preferentially  to  final  states  in  which  there  are  no 
isolated  fermions.  This  leads  to  the  selection  rule  of  BBC  reactions:  "bosons 
in  and  bosons  out".  Also,  because  the  BBC  is  composed  of  bosons,  at  low 
temperatures,  collectively  all  of  its  constituents  will  occupy  the  lowest  k=0 
energy  band,  and  a  "true"  Bose  condensate  will  form,  which  will  exhibit 
ionically  mediated  superconductivity  and  the  Meissner  effect. 

The  host  lattice  determines  the  spatial  extent  of  the  single  body  boson 
wave  function.  The  lattice  provides  periodic  potential  wells  in  which  D+  ions 
reside.  The  magnitude  of  the  zero  point  motion  rz  p  associated  with  the 
classical  turning  points  of  each  potential  well  defines  the  minimum  "zero  point 
volume"  V,  p  «  4*/3r, p 3  occupied  by  a  single  boson.  As  a  consequence,  the 
maximum  charge  density  of  each  BBC  D*  ion  is  greatly  reduced  from  its  peak  free 
space  particle  value.  We  will  show  that  the  value  of  the  lower  bound  of  c  that 
is  needed  to  produce  a  fusion  power  density  of  10  W/cm3  is  about  10'7  D+  per 
host  atom.  Since  this  value  is  much  less  than  the  D  concentration  used  in 
fusion  experiments1,  our  model  predicts  that  most  of  the  D  content  is  not  part 
of  the  BBC. 

A  key  property  of  a  BBC  is  its  ability  to  form  multiple  occupancy  band 
states.  Consider  the  idealization  of  a  group  of  bosons  which  do  not  interact 
with  each  other  but  do  interact  with  the  lattice.  These  bosons  distribute 
themselves  throughout  the  lattice  as  a  perfect  gas.  Their  distribution 
includes  the  multiple  occupancy  of  a  unit  cell.  If  the  concentration  of  these 
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bosons  is  cB,  the  concentration  of  n-boson  unit  cell  occupancies  is  cBn  in  the 
low  concentration  limit.  These  multiple  occupancies  manifest  themselves  in  the 
many  body  wave  function  as  products  of  Bloch  state  wave  functions  of  the  form 

TT’*,iBioeb»  which  may  be  viewed  as  a  multiple-particle  generalization  of  a  band 

state.  The  "4-particle  band"  is  the  probable  reactive  component  in  condensed 
matter  fusion. 

DYNAMICS  OF  BBC  SELF- INTERACTION 

In  a  BBC,  there  are  no  individual  particles.  There  are  only  collective 
states.  These  collective  states  have  properties  common  to  all  quantum  fields 
and  are  subject  to  self-fluctuation.  Because  of  the  indistinguishabil ity  of 
the  bosons  in  the  absence  of  outside  perturbations,  we  have  two  choices:  either 
1)  an  individual  boson  at  a  given  lattice  site  can  be  viewed  as  a  composite 
occupation  of  many  single  particle  eigenfunctions,  where  each  eigenfunction 
contributes  a  fraction  of  the  over-all  charge  (and  wave  function  amplitude)  and 
each  is  normalized  to  the  volume  of  the  entire  crystal,  or  2)  each  function  can 
be  normalized  to  a  single  unit  cell  whose  identity  with  respect  to  other  unit 
cells  within  the  crystal  is  lost.  The  second  option  provides  a  method  for 
identifying  individual  virtual  bosons  within  a  unit  cell.  In  either  case,  once 
a  perturbation  that  disrupts  the  periodicity  of  the  potential  is  introduced  at 
an  isolated  site,  a  proper  accounting  of  the  associated  many  body  interactions 
throughout  the  crystal  requires  that  the  initial  many  body  BBC  state  be  viewed 
in  the  framework  of  a  condensate,  made  up  of  the  products  of  1,  2,  3,  4,...,NB 
single  particle  Bloch  state  wave  functions. 

The  second  option  provides  an  alternative  to  treating  the  lattice  as  a 
whole.  It  provides  a  picture  in  which  many  "particle-like"  objects,  such  as 
di -particles  (involving  the  product  of  two  single  particle  wave  functions),  may 
participate  in  a  disruption  of  periodic  order  at  a  specific  site.  The 
"multiple-particle  bands"  that  correspond  to  these  multiple-particle-like 
objects  have  equal  density  at  each  lattice  site.  In  the  absence  of  self¬ 
interaction,  the  "multiple-particle  bands"  have  the  same  energy  as  the  sum  of 
the  individual  "single-particle  band”  energies.  Each  particle-like  object 
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(single  or  multiple)  may  be  expressed  in  terms  of  sums  of  products  of  Wannier 
states  (&/S),  where  the  density  j^J2  describes  the  transient  integer  occupancy 
at  the  n111  lattice  site.  (Each  of  these  is  a  virtual  state.)  The  virtual 
single  occupancy  state  has  a  self-interaction  energy  approximately  equal  to  the 
energy  of  a  real  D+  ion  in  the  host  lattice  since  the  BBC  is  in  thermal 
equil ibrium  with  real  interstitial  D+.  The  multiple  n-occupancy  virtual  states 
have  a  higher  electrostatic  self-interaction  energy  aE(i1)  calculable  from  their 
overlapping  charge  distributions  at  a  site.  For  the  double  occupancy  case, 
corresponding  to  the  difference  in  self-interaction  energies  between  a  doubly 
and  singly  charged  particle, 

aE<2)  -  3/2  e2  J  J  o(r)  dr')  dr  dr',  (3) 

|r-r'| 


where  />(r)  is  the  particle  density  of  the  two  particle  Wannier  state.  We  use 
a  Gaussian  to  approximate  the  ground  state  spatial  distribution  of  each  Wannier 
state  density, 

UJ2  -  (2A)3/2  exp(-2|r-RJ2/r,.p.2)/rs#Pi3  ,  (4) 

about  the  the  lattice  site  R„.  For  a  typical  value  of  rz  p  =0.5  bohr,  aEc2>  = 
81.6  eV  and  aE(4)  »  408  eV.  The  quad-deuteron  value  is  very  much  less  than  the 
self-interaction  energy  (AE(4>fr.#)  of  4  free-space  deuterons  located  at  a  common 
origin.  is  tens  of  millions  of  eV  since  the  electrostatic  self¬ 
interaction  energy  is  comparable  to  the  nuclear  self-interaction  energy  AE(4)nuc. 
Thus,  the  lattice- induced  broadening  of  the  charge  distribution  reduces  aE4  so 
strongly  that  electrostatic  effects  may  be  treated  as  a  small  perturbation 
relative  to  those  due  to  AE<4)nuc.  As  a  consequence,  it  is  possible  to  isolate 
the  effects  of  electrostatic  self-interaction  from  those  due  to  changes  in  the 
strong  potential  AV,trong  through  a  decoupling  of  the  center  of  mass  motion  from 
the  motion  associated  with  variations  in  proton-neutron  separation.  A  similar 
decomposition  is  not  possible  in  free  space  where  both  forms  of  self¬ 
interaction  proceed  at  comparable  rates.  The  electrostatic  lifetime  r(n)  of  the 
virtual  states  is  determined  by  Planck  relationship  r(n)  =  h/AE(n),  where  h  is 
Planck's  constant.  For  the  di-  and  quad-  occupancy  states,  r(2>  =  5.1  x  10'17 
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s,  and  r<4)  -  1.0  x  10'17  s,  which  should  be  compared  with  a  nuclear  Planck 
lifetime  rnuc  -  10‘22  s.  The  large  difference  in  lifetime  means  that  an 
equilibrium  condition  can  exist  between  a  spatially  compact  nuclear  state 
involving  four  deuterons  and  the  four  particle  virtual  state  over  a  timescale 
that  is  long  with  respect  to  the  nuclear  interaction  time  but  short  with 
respect  to  self-induced  electrostatic  fluctuation.  Thus,  it  becomes  possible 
for  all  of  the  ions  to  participate  collectively  in  nuclear  interactions  at 
isolated  locations. 

FUSION  FROM  THE  BBC 

Let  us  now  consider  the  fusion  problem.  The  most  probable  reaction 
appears  to  be  the  pure  bosonic,  momentum  balanced  reaction 

4D+Bloch  -  8Be++++*Bloch  -  2  Alpha ,  (5) 

or  its  Wannier  state  equivalent, 

4D  virtual  8^e  Virtual  2  Alpha  .  (6) 

In  Eqs.  (5)  and  (6),  4D+  refers  to  the  4-fold  occupancy  state  or  quad-deuteron, 
and  ®Be+*++*  stands  for  a  compact  configuration  at  nuclear  density  which 
possesses  the  symmetries  and  internal  quantum  numbers  associated  with  the 
initial  four  deuterons.  This  intermediate  state  has  the  same  total  energy  as 
the  initial  four  deuterons,  but  nucleons  within  the  "nuclear  force  bag"  have 
high  kinetic  energies.  The  final  state  requires  rearrangement  of  nucleons  so 
as  to  produce  two  spin=0  Alpha's  each  containing  23.8  MeV  relative  to  the 
center  of  mass.  The  full  reaction  is  in  principle  reversible  and  satisfies  the 
"bosons  in  and  bosons  out"  selection  rule  (all  components  can  exist  as  BBC 
states) . 

We  also  consider  the  possibility  of  two  other  reactions,  suggested  by  the 
apparent  requirement  that  LiOD  (as  opposed  to  NaOD)  be  electrolyzed  to  produce 
anomalous  heat  in  a  Pd  electrode.  With  overvoltage  electrolysis  some  Li  can 
be  expected  to  enter  the  host  lattice.  Although  the  effect  of  interstitial  Li 
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is  likely  an  electronic  structure  effect,  there  could  be  possible  nuclear 
n  cions: 


26Li%loch  -  12C~*Blocb  -  3  Alpha  +  4e’,  (7) 

and 

D+Bioch  +  6Li+i«Purity  *  aBe++*llipurit,  -  2  Alpha  +  2e*  (8) 

In  Eq.  (7),  26Li+  is  a  double  occupancy  6Li+  Bloch  state.  In  Eq.  (8),  6Li+impurity 
is  a  real  interstitial  Li+  ion  which  is  assumed  to  react  with  the  D+  BBC  initial 
state.  Equation  (8)  reflects  an  apparently  different  type  of  BBC  reaction. 
However,  since  the  final  Alpha  particle  states  can  be  considered  impurity 
states,  the  difference  may  be  more  apparent  than  real.  We  will  not  further 
consider  either  Li  reaction. 

Let  us  now  consider  whether  the  formation  of  a  D+  BBC  can  lead  to  the 
observed  heat  release  in  Pd  by  fusion  reaction  5.  Details  of  our  analysis  are 
given  in  the  Appendix  and  elsewhere3.  Using  the  observed1  value  of  the  power 
density  Pfu,ion“l°  W/cm3  (4  MJ/120  hours),  we  set  an  approximate  lower  bound 
(c1"")  for  the  concentration  (  c  )  of  injected  D+  ions,  required  to  obtain  the 
observed1  heating.  Our  estimate  of  ciow,r  follows  from  an  estimate  of  the  fusion 
rate,  rfu,ion,  based  on  the  Fermi  Golden  rule. 

Because  the  self-induced  perturbation,  AV„tron8,  associated  with  Eq.  (5) 
is  short-ranged  and  involves  interactions  between  eight  nucleons,  two 
consequences  follow  from  Eq.  (1).  1)  Large  contributions  to  the  associated 
perturbation  matrix  element  result  from  the  product  between  any  four  single 
particle  D+  Bloch  states  throughout  the  lattice  in  the  initial  state  and  the 
product  of  the  highly  localized  perturbation  with  the  final  state  functions 
associated  with  the  recoiling  alpha  particles.  2)  All  contributions  to  the 
matrix  element  are  negligible  except  for  those  associated  with  the  integrations 
over  regions  near  the  site  where  interaction  takes  place.  Then,  the  number  of 
distinct  contributions  that  arise  from  the  overlap  of  each  of  the  NB  (»1) 
single  particle  Bloch  states  associated  with  the  initial  state  (in  which  there 
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is  at  most  1  D+  wave  function  centered  at  each  site)  with  the  product  of  a Vstrong 
and  the  4  D+  -double-alpha  particle  final  state  wave  functions  near  the  location 
Rtu.ion  of  the  perturbation  is  Nb(Nb-l)(Nb-2)(Nb-3)  -  Nb\  Hence,  given  an 
approximate  value  for  the  energy  release  per  fusion,  Qfusion,  for  fixed  values 
PfU*ion  *  ^"fusion  Qfusion>  ^ t  follows  that  Pfuli0n  a  £  ^nue  »  where  rnuc  is  the 
reaction  rate  of  a  highly  localized  nuclear  reaction  at  the  single  site  Rf,.siotl. 

To  evaluate  clow,r  we  determine  an  upper  bound  for  rnuc  using  a  simplified 
reaction  model  for  Eq.  (5).  This  simplified  model  provides  a  higher  reaction 
rate  than  is  expected.  We  replace  the  eight  nucleon  problem  of  the  reaction 
by  an  independent  particle  nucleon  model*  in  which  the  fusion  process  is 
modelled  through  a  square  well  perturbation  of  the  nucleonic  potential.  This 
square  well  is  defined  by  a  suitable  upward  shift  in  the  kinetic  energy  T0  of 
each  D+  ion  over  a  characteristic  nuclear  dimension  within  the  vicinity  of 
Rfusion *  This  is  equivalent  to  requiring  in  the  reverse  of  reaction  5,  an 
immediate  dissipation  of  the  energy  of  both  alpha  particles  (which  may  be 
viewed  as  a  short-lived  virtual  state)  in  the  immediate  region  in  which  the 
nuclear  reaction  takes  place.  The  height  of  the  square  well  is  Qfusion  =  47.6 
MeV  (calculated  from  the  decrease  in  rest  mass5),  which  is  inferred  from  the 
assumption  that  the  wave  function  overlap  associated  with  the  quadruple 
deuteron  "multiple  particle  like  object"  may  be  viewed  as  an  excited  state  of 
a  8Be++++  nucleus. 

Explicitly,  using  a  suitable  set  of  minimal  uncertainty  gaussian  wave 
packets  in  which  the  initial  state  is  defined  with  a  characteristic 
electrostatic  dimension  rt  p  =  0.26  A  and  volume  Vz  p  =  4/3  *  rt  p  \  while  the 
final  state  is  confined  to  the  nuclear  volume,  Vnuc  =  4/3  *-  RB.3,  defined  by  the 
Be****  nuclear  radius  R„,  *  9.108  x  10'13  cm,  we  find  that 


T*.im  ■  [2  «d  R».V(3H2)]  (16/3)3'!  [Vnuc/Vj  p  ]QfuBlon2 


X 


(clow,r)4  , 


(9) 
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where  Mo  is  the  mass  of  a  deuteron  and  K=h/27r.  Using  P£uilon  =  10  W/cm3  =  Nslt. 
Q fusion  r£u.io,,  where  Nllt.  (s  4.81  x  1022  cm'3)  is  the  density  of  Pd,  and  Eq.  (9), 
we  find  that  £low,r  equals  2.8  x  10'7.  Though  Eq.  (9)  was  derived  for  an 
isolated  perturbation,  which  applies  to  low  temperature,  the  identical  argument 
can  be  gener'.ized  to  the  case  of  higher  temperature  by  performing  a  suitable 
averaging  over  initial  state  band  energies  in  the  Fermi  Golden  rule.  In  this 
case,  it  is  still  true  that  P£usion  «  c4  rnuc.  However,  rnuc  and  c*  may  both 
become  strongly  dependent  on  temperature  due  to  disruption  of  Bloch  symmetry 
by  the  fusion  process.  This  feature  could  prove  valuable  for  controlling  the 
process. 

MANAGING  FUSION  THROUGH  THE  BBC 

The  BBC  D+  fusion  picture  requires  a  suitable  host  lattice.  One  critical 
feature  enabling  condensed  matter  fusion  is  the  establishment  of  conditions  in 
which  some  D  resides  in  the  lattice  in  ionic  form  screened  by  itinerant 
electrons.  A  second  requirement  is  that  energies  of  the  reactant  single 
particle  band  be  close  to  that  of  chemically  bonded  D  so  that  it  is  possible 
to  populate  BBC  states  at  realizable  temperature.  A  third  requirement  is  that 
the  lattice  provide  adequate  regularity  and  size  to  maintain  approximate  Bloch 
symmetry.  A  large  rz  p  appears  favorable  and  reflects  a  low  interstitial  D+ 
binding  energy,  important  for  reducing  electrostatic  self-interaction.  A  large 
r*.P.  favors  high  diffusivity.  It  is  interesting  that  Pd  provides  a  diffusivity6 
for  D  that  is  two  times  greater  than  for  H,  contrary  to  expectation,  suggesting 
that  the  bosonic  character  D*  may  contribute  substantially  to  the  diffusivity. 
This  line  of  reasoning  suggests  that  a  Pd-Ag  host,  which  provides  enhanced  H- 
diffusion,  could  be  favorable  for  fusion  reactions. 

The  key  observable  of  BBC  fusion  is  intense  Alpha  emission.  For  the  4- 
deuteron  reaction  (Eq . (5) ) ,  the  Alpha  particles  have  a  range  in  Ilford  emulsion 
of  more  than  240  microns5.  Strong  M-  and  N-  shell  x-ray  line  emission  at  -0.68 
and  -0.094  keV  should  be  present.  Neutrons  are  not  predicted  except  as  a 
result  of  Alpha  bombardment  of  host  impurities  or  of  the  surrounding 
electrolysis  medium.  Ionic  superconductivity  at  sufficiently  low  temperature 
is  expected. 
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APPENDIX 


DERIVATION  OF  EQUATION  (9) 

In  this  appendix,  we  provide  additional  detail  concerning  the  derivation 
of  Eq.  (9).  We  evaluate  r£usion  from  the  Fermi  Golden  rule: 

rfu.ion  »  2V  h  s£  s(E£  -  Ei)|<*£|AV|*iBl0Ch>|2  .  (A-l) 


Here,  the  many  body  Bloch  bosonic  initial  state  involves  a  suitably 
symmetrized  sum  of  products  of  single  particle  wave  functions.  Each  product 
consists  of  a  macroscopic  number,  Nb,  of  occupied  single  particle  Bloch  wave 
functions  associated  with  a  particular  distribution  of  Nb  distinguishable  D+ 
ions,  in  which  not  more  than  one  single  particle  bosonic  Bloch  wave  function 
^Bioch  centere(j  a  specific  site,  and  each  function  has  an  associated 
energy  band  eigenvalue,  from  the  set  {et,  i«=l,Nb)  of  occupied  bands.  The  sum 
over  products  is  used  to  impose  Bose  symmetry,  so  that  *iBl0Ch  is  invariant  when 
the  coordinates  of  any  number  of  its  single  particle  wave  functions  are 
interchanged.  Since  each  ion  is  located  at  one  of  Nb  different  sites,  c  = 
Nb/NL>tt»  where  NL#tt  is  the  total  number  of  sites  in  the  crystal .  For  the 
relevant  case  in  which  £  is  small,  each  D+  wave  function  ^iBloch  is  centered  at 
a  lattice  site  far  from  the  other  lattice  sites. 

The  many  body  final  state  wave  function  »£  of  Eq.  (A-l)  consists  of  a 
similar  sum  of  products  of  single  particle  wave  functions,  but  for  4  fewer 
single  particle  Bloch  states.  In  this  case,  Nb-4  band  energies  and  wave  vectors 
and  the  associated  wave  functions  are  selected,  once,  and  for  all,  to  be 
identical  to  Nb-4  eigenvalues  from  the  initial  state.  Each  term  in  the  final 
state  sum  consists  of  a  product  of  the  associated  Nb-4  initial  state  single 
particle  wave  functions,  centered  at  different  sites,  and  4  single  "defect" 
D+  wave  functions,  which  are  all  centered  at  the  site  Rfusion.  Also  in  Eq.  (A- 
1),  h  is  h/(2*),  aV  is  the  self-interaction  perturbation  of  the  periodic 
potential . 
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Because  of  the  important  lattice-induced  broadening  of  the  charge  and 
reduction  in  the  electrostatic  self-interaction  energies  AE<n)  for  the  virtual 
occupation  by  an  n-particle  Wannier  state  fluctuation,  we  may  approximate  aV 
by  AV,trons.  At  low  temperature  (T),  only  isolated  nucleonic  perturbations  need 
be  considered  because  adjacent  or  near-by  perturbations  result  in  increases  in 
entropy  and  free  energy  of  the  crystal  system.  Near  T=0,  there  is  only  a 
single  site  (Rfu,ion)  where  AV,trons  is  non-vanishing.  Then,  in  the  evaluation  of 
the  transition  matrix  element  <MAV»tr0ngl'i,iBl0Ch>»  because  AV8trong  is  short- 
ranged  and  each  single  particle  initial  state  wave  function  possesses  Bloch 
symmetry,  two  important  conseqences  follow  from  Eq.  (1).  1)  Large  contributions 
to  the  perturbation  matrix  element  result  from  the  product  consisting  of  any 
four  single  particle  Bloch  states  throughout  the  lattice  in  the  initial  state 
multiplied  by  the  highly  localized  perturbation  with  the  four  defect  final 
state  i>t  functions.  2)  All  contributions  to  the  matrix  element  are  negligible 
except  for  those  associated  with  the  integrations  over  regions  near  the  site 
where  interaction  takes  place. 


The  number  of  distinct  contributions  that  arise  from  the  overlap  of  each 
of  the  Bloch  states  with  the  4  defect  states  in  the  region  near  Rfusion  is  Nb(Nb- 
1 ) (Nb-2) (Nb-3)  -  Nb4,  for  Nb  »  1.  Then,  when  each  single  particle  Bloch  state 
is  normalized  to  the  volume  of  the  crystal,  it  follows  as  a  consequence  of  Eqs. 
(1)  and  (A- 1 )  that 

2,/  ft  tZ,.(|4V“”VJP,t(eJ)  )]  £*  •  (A-2) 

Here  AVstton*£.  j  is  the  matrix  element  derived  from  integration  over  the  region 
of  volume  Vnuc,  centered  at  R£uIlon,  where  the  product  of  a V8trong  with  the  four 
localized  final  state  single  particle  wave  functions  and  four  initial  state 
functions  is  non-negl igible.  Here,  p£(e j)  is  the  final  density  of  states 
evaluated  with  energy  ej-Sj  ej,  where  the  summation  includes  the  four  initial 
state  bands  that  are  excluded  from  the  final  state. 

As  discussed  in  the  text,  to  establish  cl0",r  in  Eq.  (9),  it  is  sufficient 
to  consider  a  faster,  simpler  reaction  based  on  Eq.  (5):  we  have  replaced  the 
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eight  nucleon  problem  of  the  reaction  by  an  independent  particle  nucleon  model4 
in  which  the  fusion  process  is  modelled  through  a  square  well  perturbation  of 
the  nucleonic  potential.  This  square  well  is  defined  by  a  suitable  upward 
shift  in  the  kinetic  energy  Te  of  each  D+  ion  over  a  characteristic  nuclear 
dimension  within  the  vicinity  of  Rtullon.  This  is  equivalent  to  requiring 
immediate  dissipation  of  the  energy  of  both  alpha  particles  (which  may  be 
viewed  as  a  short-lived  virtual  state)  in  the  reversed  reaction  in  the 
immediate  region  in  which  the  nuclear  reaction  takes  place. 

Then,  each  defect  and  initial  state  function  eigenvalue  is  measured 
relative  to  the  constant  energy  scale  associated  with  a  perfectly  ordered 
environment  in  which  no  fusion  occurs.  The  height  of  the  square  well  is  Qfu,ion 
-  47.6  MeV,  which  is  inferred  from  the  assumption  that  the  wave  function 
overlap  associated  with  the  quadruple  deuteron  "multiple  particle  like  object" 
may  be  viewed  as  an  excited  state  of  a  ®Be++++  nucleus  and  from  the  associated 
decrease  in  rest  mass  due  to  the  decay  of  the  unstable  ®Be++++  nucleus  to  two 
alpha  particles. 

Explicitly,  we  find 

rfu.ion  -  2*/X  (  |0{J)j2  pt(ej)  )  c*  Qfu,ion2  ,  (A-3) 

where  the  overlap  matrix  element  0(J)  is  given  by 

0(J)  -  TTr«,d  j  -  i.*  J„u=  *r‘(r')  *jBloch(r')  dV,  (A-4) 

and  each  integral  in  Eq.  (A-4)  extends  only  over  the  region  near  RCulion  where 
the  values  of  AV,tton8  associated  with  the  square  well  are  non-negl igible. 
Because  AV.1>e  «  AV,trons,  we  may  rewrite  each  product  of  four  initial  and  final 
state  wave  functions  in  a  separable,  Born-Oppenheimer  form,  involving  a  slowly 
varying  envelope  function  fa  whose  spatial  dependence  only  involves  the 
center  of  mass  coordinate  R,,,,  {famfa(Rca)) ,  multiplied  by  three  more  rapidly 
varying  functions,  ifrnue,  whose  coordinate  dependencies  are  all  expressed 
relative  to  R,,,.  Then,  we  obtain  a  reasonable  upper  bound  for  |0(J)|2  pf(ej)  by 
considering  appropriate  bounds  for  the  single  integral  over  the  center  of  mass 
coordinate  of  the  product  fait  fa  L  and  including  only  the  dependence  of  the 
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associated  change  in  total  energy  in  pf(e j)  using  the  single  particle  density 
of  states  that  is  appropriate  for  VtM.f  Thus,  we  find 


0(j)  -  Jnuc  **.<>')  ^cw.iBloeh(^)  d3r'  ,  (A-5) 

where  the  single  integration  extends  over  the  volume  Vnue  centered  at  Rfusion 
because  AV,tron<  is  non-vanishing  only  in  this  region.  Then,  the  maximum  in  0(J) 
occurs  when  ^  f  is  primarily  localized  within  Vnuc  but  has  considerable  overlap 
with  This  occurs  when  f  is  appreciable  only  on  length  scales 

associated  with  nuclear  interaction  while  Vcm.i  has  variation  on  a  considerably 
larger  length  scale. 

We  derive  £l0"#r  from  Eqs.  (A-3)-(A-5)  using  initial  and  final  states, 
approximated  from  minimal  uncertainty  wave-packets,  ifrpaek,tii  and  tfpackatjf  that 
Incorporate  this  interplay  between  length  scale.  Each  packet  possesses  a 
suitable  characteristic  variance  (<r2ehaP,et,ri,tic>)  in  position  associated  with 
the  minimal  error  in  the  measurement  of  position  and  momentum.  Specifically, 

tfp.ck.t.1  -  Nc  exp(-r2/r,  p  2)  ,  (A-6) 

where  N„  -  (2/(irr,  p 2))3/\  and 

^paek«t . t  -  N f  exp(-3r2/(4  Rb.2))  ,  (A-7) 

and  Nf«  (3/(2wRBa2))3/4.  Here,  RBa*  9.108  x  10"13  cm,  as  derived  from  a  second 
calculation  associated  with  the  independent  nucleon  model  of  Be  in  which  all 
of  the  nucleons  are  bound  by  a  square  well  of  depth  Vo*50  MeV,  the  number  of 
bound  nucleon  states  is  precisely  9,  and  the  associated  Fermi  energy  for  the 
eight  nucleon  case  (-17  MeV  below  the  top  of  the  well)  is  typical  for  nuclear 
matter4.  Also,  -  4/3wRBa3.  Since  in  Eq.  (A-6),  realistic  values  of  rt  p  can 
be  as  large  as  a  significant  percentage  of  the  nearest  neighbor  Pd  separation 
Aw  (2.75  A),  we  have  used  rt  p  -  0.5  B  «  0.265  A  *  0.1  APd.  Then,  Vnuc/VI  p 
-  10"“,  where  V,  p  -  4/3  wra  p  3  is  the  total  volume  of  each  D+  due  to  zero  point 
motion.  From  Eqs.  (A-5) -(A-7),  we  find 
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0(J)  -  (16/3)3'*  (Vnu=  /  )*« 


i 


(A-8) 


where  we  have  approximated  4RB.2+3rI  p  2  with  3rz  p  2  since  RB,2  «  r2  p  2.  The 
appropriate  final  density  of  states  is 

Pt( 0)  -  1/(H  «nuc)  ,  (A-9) 

where  -  3/2  ft  /  (Mp  RB#2).  Equation  (A-9)  follows  from  Eq.  (A-7)  because 

^p«ck#t.f  the  ground  state  of  a  three  dimensional  harmonic  oscillator  that 

possesses  an  angular  frequency  of  wnue.  Substituting  Eqs.  (A-8)  and  (A-9)  into 
Eq.  (A-l),  we  find 

IW™  -  2/3  .  (16/3)3'2  {Qfull„3  („  R,.  )3 

X  V„UI)/(2  H  V,.t.  V„)  (s1”")*,  (A-9) 

where  *  -  (Vn  Mo  2/fi2)1'2 ,  which  is  equivalent  to  Eq.  (9). 
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